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ABSTRACT 

We present deep, 1.6/rm, continuum images of eleven high-redshift (0.811 < z < 1.875) 3CR radio 
galaxies observed with NICMOS on-board the Hubble Space Telescope. Our NIGMOS images probe 
the rest-frame optical light where stars are expected to dominate the galaxy luminosity. The rest-frame 
ultraviolet light of eight of these galaxies demonstrates the well-known “alignment effect”, with extended 
and often complex morphologies elongated along an axis close to that of the FRH radio source. As has 
been previously noted from ground-based near-infrared imaging, most of the radio galaxies have rounder, 
more symmetric morphologies at rest-frame optical wavelengths. Here we show the most direct evidence 
that in most cases the stellar hosts are normal elliptical galaxies with law light profiles. For a few 
galaxies very faint traces (less than 4% of the total 7J-band light) of the UV-bright aligned component 
are also visible in the infrared images. We derive both the effective radius and surface brightness for nine 
of eleven sample galaxies by fitting one- and two-dimensional surface-brightness models to them. We 
compare the high-redshift radio galaxies to lower redshift counterparts. We find their sizes are similar 
to those of local FRH radio source hosts and are in general larger than other local galaxies. The derived 
host galaxy luminosities are very high and lie at the bright end of luminosity functions constructed at 
similar redshifts. This indicates that the high-redshift radio galaxies are likely rare, massive sources. The 
galaxies in our sample are also brighter than the rest-frame size-surface-brightness locus defined by the 
low-redshift sources. Passive evolution roughly aligns the z w 1 galaxies with the low-redshift samples 
with a slope equal to 4.7. This value falls intermediate between the canonical Kormendy relation (pe 3.5) 
and a constant luminosity line (= 5). The optical host is sometimes centered on a local minimum in 
the rest-frame UV emission, suggesting the presence of substantial dust obscuration. We also see good 
evidence of nuclear point sources (no brighter than 5% of the total 7J-band light) in three galaxies. 
Overall, our results are consistent with the hypothesis that these galaxies have already formed the bulk 
of their stars at redshifts greater than z > 2, and that the AGN phenomenon takes place within otherwise 
normal, perhaps passively evolving, galaxies. 

Subject headings: galaxies: formation — galaxies: active — infrared: galaxies 


1. INTRODUCTION 

A variety of techniques have been used to find galaxies 
at high redshift, but the hosts of powerful radio sources 
remain uniquely interesting as the most luminous and ar¬ 
guably the most massive galaxies known at z > 1. Ra¬ 
dio surveys have successfully found many distant galaxies 
with bright optical or infrared counterparts which, in addi¬ 
tion, often have prominent emission-lines. As the number 
of known radio galaxies (RGs) increased, it was discov¬ 
ered that selection via radio flux was preferentially find¬ 
ing early-type galaxies. At low to intermediate redshift, 
powerful radio sources are almost exclusively found in gi- 

^ Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astron¬ 
omy, Inc., under NASA contract NAS 5-26555. These observations 
are associated with proposal #7454. 

^ National Optical Astronomy Observatory is operated by the As¬ 
sociation of Universities for Research in Astronomy, Inc. (AURA), 
under cooperative agreement with the National Science Foundation. 


ant elliptical (gE) host galaxies, which are often situated 
in moderately rich groups or clusters (Matthews, Mor¬ 
gan, & Schmidt 1964; Hill & Lilly 1991). Several radio 
galaxies at high-redshift have also been found to inhabit 
overdense environments; a sign that RGs mark deep po¬ 
tential wells. These facts suggest that high-redshift radio 
galaxies (HzRGs) are the progenitors of present-day gi¬ 
ant galaxies. In addition, the recent discovery of strong 
correlations between galaxy and central black hole masses 
(Magorrian et al. 1998; Gebhardt et al. 2000; Ferrarese & 
Merritt 2000) implies an evolutionary scenario in which 
all giant galaxies at all observed redshifts harbor super- 
massive black-holes which undergo (perhaps recurrent) ac¬ 
cretion episodes to the present day. If not prone to strong 
selection effects, these relatively rare sources make excel¬ 
lent discriminators of structure formation and evolution 
models. However, the question remains: how might the 
presence of an AGN (both dynamically and its associated 
radiation and outflows) affect a host galaxy’s formation 


1 


2 


Massive Ellipticals at z w 1 


and evolution, particularly in the extreme cases such as 
the SCR radio sources? 

There is evidence that the powerful nuclei of radio galax¬ 
ies play a major role in the stellar hosts’ evolution. Local 
radio galaxies frequently show morphological peculiarities 
such as tidal tails and dust lanes (e.g., Heckman et al. 
1986; Smith & Heckman 1989) which are usually absent 
in the general galaxy population. While these features do 
not usually dominate the galaxy light (i.e., they can still 
be classified as ellipticals), they have been interpreted by 
some to suggest that radio activity is triggered by galaxy 
interactions and mergers (e.g., Heckman, Carty, & Bothun 
1985). At higher redshifts, z > 0.7, powerful radio galax¬ 
ies have strong rest-frame UV line and continuum emission 
which is often extended on an axis similar to that of the 
radio source (McCarthy et al. 1987; Chambers, Miley, & 
van Breugel 1987). This “alignment effect” has a domi¬ 
nant influence on the measured morphological and photo¬ 
metric properties of high-redshift radio galaxies, at least 
at ultraviolet wavelengths where the blue, AGN-related 
continuum emission outshines the stellar light from the 
presumed extant host galaxy. 

Optical observations with H5T/WFPC2 have revealed 
the remarkably peculiar and complex morphologies of this 
aligned light (Dickinson, Dey, & Spinrad 1995; Best, Lon- 
gair, & Rottgering 1997; McCarthy et al. 1997). The 
rest-frame UV continuum is often segregated into several 
discrete clumps which are strung along the radio axis. 
There are three viable explanations for the alignment ef¬ 
fect: scattering of AGN light by dust or electrons, star- 
formation induced by the passage of the radio jets through 
the galactic ISM, nebular continuum emission (Dickson 
et al. 1995) or some combination thereof. It has been es¬ 
tablished from the high degree of polarization that the spa¬ 
tially extended, aligned, UV continuum emission in many 
2 « 1 — 2 RGs is scattered light (e.g., Jannuzi & Elston 
1991; di Serego Alighieri, Gimatti, & Fosbury 1994; Dey 
et al. 1996; Gimatti et al. 1997; Tran et al. 1998). Data on 
the higher redshift RGs {z > 2) are more limited, but at 
least one galaxy shows evidence of young stars contribut¬ 
ing to the aligned light (Dey et al. 1997). If alignment by 
induced star-formation is prevalent at early times, then a 
large fraction of the stars in the galaxy may have been 
formed via this process. 

On the other hand, near-infrared observations of distant 
radio galaxies have found a smooth magnitude-redshift re¬ 
lation {K-z) with remarkably small scatter out to at least 
2 w 3 (e.g., Lilly & Longair 1984; Eales et al. 1993; Mc- 
Garthy 1993; Jarvis et al. 2001). The relation is well- 
described by a coeval, passively evolving galaxy popula¬ 
tion which formed at z > 4. Less luminous radio galaxies 
(Eales et al. 1997) over the same redshift range display 
a K-z relation analogous to their powerful SCR counter¬ 
parts. In effect the radio selection has chosen a particular 
class of galaxy from the general population. This is consis¬ 
tent with the local association of radio sources with giant 
elliptical galaxies. 

To explain the alignment effect and K-z diagram simul¬ 
taneously, previous authors (Rigler et al. 1992; Best, Lon¬ 
gair, & Rottgering 1998) have adopted a two-component 
source model. A massive, old, early-type galaxy accounts 
for the K-z relation and a bright, but dynamically insignif¬ 


icant, AGN-related contribution (scattered nuclear light 
and/or jet-induced star formation) dominates the rest- 
frame UV light. Ground-based infrared observations have 
supported this view (Rigler et al. 1992; Best, Longair, & 
Rottgering 1997), but have been limited by the relatively 
poor angular resolution attainable from the ground. How¬ 
ever, two-dimensional fitting of higher angular resolution 
HS'T/WFPG2 imaging of several z « 1 SCR radio galaxies 
(McLure & Dunlop 2000) find elliptical-like host galaxies. 
Such direct observations of the elliptical component can 
corroborate this theory and perhaps provide a quantify 
the biases inherent in using AGN hosts to trace galaxy 
evolution. 

To target the stellar host and its morphology requires 
three observational capabilities: 1) access to rest-frame 
optical or longer wavelengths to increase the contrast be¬ 
tween stellar light and AGN-related emission, 2) filters 
which target continuum emission and avoid bright emission¬ 
lines, and 3) high angular-resolution to discern morpholog¬ 
ical features with clarity. HST/NICMOS currently pro¬ 
vides the best technology to address these three criteria. 
We therefore obtained HST/NICMOS imaging of a sample 
of 11 powerful radio galaxies at 0.8 < z < 1.8. Our targets 
are a representative subsample of the most powerful radio 
galaxies at these redshifts. Although these 11 galaxies do 
not constitute a statistically complete sample of high red¬ 
shift radio galaxies, they do span the full range of optical 
and near-infrared (NIR) photometric properties shown by 
this population as a whole. For example, the strength of 
the alignment effect in our sample ranges from essentially 
zero (e.g., 3CR 65) to very well aligned (e.g., 3GR 368 or 
3GR 266). 

We present both H5T/NICMOS and HST/WFFC2 imag¬ 
ing from our own proposals as well as from the HST data 
archive. In the next section we describe the data acqui¬ 
sition and reduction. In §3 we describe our profile fit¬ 
ting procedures and error analysis. In §4 we present the 
derived physical parameters and compare them to local 
galaxy samples, and discuss the implications of our results. 

Throughout this paper we assume Hq = 65 km s“^Mpc“^ 
and (HajIIm) = (0.7,0.3). With this cosmology 1" sub¬ 
tends 8.6 kpc at z = 1. 

2. OBSERVATIONS AND DATA REDUCTION 

The sample galaxies were chosen from the 3GR radio 
survey to represent powerful radio galaxies at z « 1. Our 
11 galaxies cover 0.8 < z < 1.8 and exhibit the range 
of host galaxy properties present in the 3GR as a whole. 
As with any flux-limited sample, 3CR galaxies show a 
strong radio-luminosity-redshift correlation. This should 
not greatly affect most of our conclusions, but should be 
kept in mind when comparing radio galaxies over different 
epochs, as we do in §4. Table A lists the sample galaxies 
(Gol. 1), their coordinates (Col. 2-5) and redshifts (Col. 6) 
along with the exposure parameters, including exposure 
time (Col. 9), filter (Col. 7), number of dithers (Col. 10), 
rest-frame wavelengths covered (Col. 8), and estimated 
emission-line contamination (Col. 11). 
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2.1. Data Acquisition 

Our NICMOS observations were taken with Camera 2 
between December 1997 and July 1998. The NICMOS 
filters F160W (1.60/rm) and F165M (1.65/rm) were used 
to obtain images in rest-frame wavelengths ranging from 
approximately 6000 to 9000A (exact wavelength ranges 
can be found in Table A). The alignment effect is es¬ 
pecially prominent in the narrow emission line gas (e.g., 
McCarthy, Spinrad, & van Breugel 1995). Therefore, to 
avoid this contamination of the stellar component in terms 
of both luminosity and morphology, we used the known 
redshifts of the sources and the composite RG spectrum 
of McCarthy (1993) to choose the filter with the least line 
contamination. We attempted to avoid potentially bright 
emission lines such as Ha and, for the higher-redshift ob¬ 
jects, [OIII] AA 4959,5007. In the case of 3CR 256 we took 
exposures in both F160W and F165M to obtain contin¬ 
uum and [OIII] AA 4959,5007 images. The resulting images 
show some contamination which has been estimated by us¬ 
ing a template spectrum with and without emission lines. 
We note that these estimates assume a spatially-constant 
equivalent width, whereas this value is likely to vary over 
the area of the galaxy. The estimated emission-line con¬ 
tamination is listed in column 11 of Table A. The tele¬ 
scope was dithered between exposures for all the objects, 
with 8-13 dither positions per object, to better sample 
the instrumental PSF, which is already critically sampled 
at 1.6/xm, and to aid with the recognition and elimina¬ 
tion of data artifacts such as dead pixels or “grot”, hot 
pixels, the unstable central column of NIC 2, and persis¬ 
tent cosmic-ray afterglow. The exposures themselves were 
taken using a combination of the MIF512 and MIF1024 
NICMOS MULTIACCUM readout sequences, providing 
exposure times of 512 and 1024 seconds, respectively. In 
Figure 1 we present the NICMOS images (second pan¬ 
els), the rest-frame UV images from WFPC2 (first panels), 
which were registered and interpolated to the NICMOS 
images, and the two-dimensional fit residuals (third pan¬ 
els), which were created by subtracting the best-fit model 
galaxy from the data as described in §3.1.2. 

2.2. Data Reduction 

The individual frames were reduced using a combina¬ 
tion of standard IRAF^/STSDAS pipeline routines (CAL- 
NICA) and our own techniques. The most prominent 
anomalies in our data were: the “pedestal”, a random 
DC offset which varies by detector quadrant; cosmic-ray 
persistence, time-decaying charge deposits from cosmic- 
ray hits during HST’s passage through the South Atlantic 
Anomaly; and noise associated with the synthetic dark 
frames in the standard pipeline (Dickinson 1999). The 
“pedestal” subtraction was performed using the PEDSKY 
task now included in STSDAS, but originally written for 
use with this dataset. It performs an iterative subtrac¬ 
tion of both a constant offset per detector quadrant and a 
“sky” component which is modeled as a constant (over the 
whole detector) multiplied by the flat-held pattern. The 
solution minimizes the pixel-to-pixel RMS in each quad- 

® IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 


rant to converge on the best fits to these parameters. 

The synthetic darks now used in CALNICA are “super¬ 
darks” generated using many more individual frames than 
the darks which were available when we reduced these 
data. Also, it is now possible to generate temperature- 
dependent dark frames for individual exposures, which 
help reduce bias gradients in the images (called “shad¬ 
ing” ). To correct our data for the noise pattern introduced 
by the original darks, we constructed median frames from 
the sky and pedestal subtracted images for each filter and 
readout sequence combination. We combined both our 
data and those from a similar observing program (M. Gi- 
avalisco, private communication) using a median, and sub¬ 
tracted the resulting image from our data. We were able 
to reduce the amount of pixel-to-pixel noise introduced by 
the artificial darks by several percent with this method. 

Fortunately, cosmic-ray persistence affected only a small 
minority of our images. The persistent pixels have charge 
deposited in them from a cosmic-ray hit, which depletes 
roughly exponentially with time. Occurrence of cosmic- 
ray persistence is strongly correlated with the HST’s pas¬ 
sage through the South Atlantic magnetic anomaly (SAA). 
Data for one of our targets, 3CR 65, were severely affected 
by post-SAA persistence, and we have used an experimen¬ 
tal method to attempt to correct for this effect. To do 
so, we “sacrificed” some of the exposures closest in time 
to the SAA passage, and used these as a model of the 
persistence pattern. We then scaled and subtracted this 
persistence model, minimizing the background RMS in a 
fashion similar to that used by PEDSKY. This worked 
at the w 10% level providing significantly improved noise 
and cosmetic characteristics. A more robust method us¬ 
ing darks automatically scheduled after SAA passages has 
been implemented for observations with NICMOS in HST 
Cycle 11. 

The individual exposures for each galaxy were combined 
using the DRIZZLE task (Fruchter & Hook 2002), with rel¬ 
ative shifts determined using image cross-correlation. In¬ 
put weight maps (inverse per-pixel variance) were gener¬ 
ated with a noise model which includes the contributions 
from the sky, readout amplifier glow, standard dark cur¬ 
rent and read noise. The DRIZZLE task then generates 
an output weight map based on the combination of input 
error maps, transforming the input in exact analogy with 
the data themselves. These output weight maps were used 
to calculate our error estimate in the subsequent analysis. 
DRIZZLE does not account for the correlation of the noise 
between adjacent pixels when it creates the output weight 
map. To correct the variance for inter-pixel correlation, 
we applied a single multiplicative factor (= 1.84) based 
on the equation derived by Casertano et al. (2000). The 
effective plate scale after drizzling is one-half the original 
scale, or 01'0375 / pixel (the diffraction-limit of HST at 
1.6/rm is « 0.168"). Several original tasks were created 
for these reductions are now included in the STSDAS re¬ 
duction package. 

The WFPC2 frames included in the present analysis are 
from CTO proposals (1070 and 6235; PI M. Longair) from 
the HST archive and several GO proposals (5429, PI W. 
van Breugel; 5925, PI P. Eisenhardt; 5967, PI M. Dickin¬ 
son; 5995, PI O. Le Fevre). Reduction of the WFPG2 im¬ 
ages was done using the standard STSDAS pipeline. The 
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WFPC2 thumbnail images presented in this work have 
been registered with the NICMOS images and interpo¬ 
lated to the scale of the drizzled NICMOS images. Ta¬ 
ble A presents magnitudes derived from the reduced data 
using aperture (Col. 3-4), total (Col. 5), and point source 
(Col. 6) photometry. These numbers have relatively little 
dependence on the model fits described in the next sec¬ 
tions. 

3. ANALYSIS 
3.1. Profile Fitting 

To quantitatively compare the morphologies of our galax¬ 
ies to other samples, we fit analytic surface-brightness mod¬ 
els to the data, primarily concentrating on exponential 
disk and -law (de Vaucouleurs 1948) profiles. The fits 
were used to determine characteristic sizes and to differ¬ 
entiate between the spheroidal and disk models. In ad¬ 
dition, a Sersic profile (r^/”-law; Sersic 1968) was also 
fit to discern any deviations from the and disk mod¬ 
els. We decided not to fit more complex models (such as 
the “Nuker” law of Lauer et al. 1995) with more free pa¬ 
rameters, because our data are unable to constrain many 
of the added parameters. For instance, the “break” ra¬ 
dius in the “Nuker” law is empirically a 100 — 300 parsecs, 
which is smaller than the physical resolution of our data. 
We used both one-dimensional (ID) and two-dimensional 
(2D) fitting algorithms to derive the best-fit profile pa¬ 


rameters. For the ID fits, the initial source position angle 
(PA) and ellipticity were fixed to measurements using the 
STSDAS task ELLIPSE. The degree of control provided 
by ELLIPSE allows easy quantification of the ellipticity 
and PA. For the 2D method, we performed fits allowing 
the PA and ellipticity to vary and also fixing them to the 
values determined from our ELLIPSE measurements. We 
subsequently used the fit with the lower reduced chi-square 
iXu) value. As it turned out, all galaxies were better fit 
using the fixed ellipticity and PA model. 

The cases of 3CR 470 {z = 1.653) and 3CR 256 {z = 
1.819) are complicated and warrant separate analyses. Our 
3CR 470 images are contaminated by the presence of strong 
Ha emission in the bandpass. It is also not clear which 
of the three sources we detect is the radio source coun¬ 
terpart. 3CR 256 is certainly the most peculiar, as well 
as the highest redshift, object in this sample. Its peak 
surface-brightness and total magnitude are much fainter 
than those of our other sample galaxies (see Table A). 
3CR 256 is also known to lie well-below the K-z relation 
formed by other powerful radio galaxies (Chokshi & Eisen- 
hardt 1991; Simpson et al. 1999). Although we present the 
imaging data, we have chosen not to include 3CR 470 nor 
3CR 256 in the current analysis. 

3.1.1. One-dimensional Fits 

To fit analytic profiles to the data, we first constructed 
azimuthally-averaged profiles (i.e., one-dimensional func- 
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Fig. 1.— WFPC2, NICMOS, and model residual images for our 9 fit galaxies, presented in ascending redshift order (from SCR 265 at 
z = 0.811 to SCR 256 at z = 1.819). The images are shown in the standard North up and East left orientation. Note the similarity between the 
NICMOS residuals and the rest-frame UV WFPC2 images, suggesting a simple connection between the two, i.e., that the NICMOS residuals 
are simply the long-wavelength tail of the light seen in the UV. The data for SCR 470 and 256 are presented despite not being analyzed in 
the same manner. The contours are at 5, 10, 20 and lOOcr above (or below) the locally determined background. The major tickmarks are 0'.'5 
apart and the panels are S7S75 on a side 
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tions of semi-major axis) for both the real data and a 
set of artificial galaxies. The analytic galaxy templates 
span a wide range of effective radii (a grid of 200 bins 
ranging from 10 — 150 drizzled pixels, or 0.375 — 5.625") 
and ellipticity (a grid of 5 bins ranging from 0.05 - 0.7). 
These were convolved in two dimensions with a simulated 
NICMOS Camera 2 PSF. We used the TinyTim program 
(Krist 1995) to generate this PSF on a sub-sampled grid 
corresponding to our final DRIZZLE pixel scale. 



0 0.5 1 1.5 
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Fig. 2. — Azimuthally-averaged intensity versus radius for both 
the TinyTim model PSF (circles) and a bright star in the field of 
SCR 470 (stars). The two curves match one another well at r < 0f^5. 
The discrepancy at larger radii will only affect the photometry of 
bright point sources, i.e., ones much brighter than those found in 
our sample. 

We then smoothed this PSF with a boxcar filter rep¬ 
resenting the original pixel scale to emulate the effects of 
under-sampling due to the finite pixel size. The resulting 
PSF has approximately the same flux distribution as the 
final reduced data. Optimally, we would have liked to use 
an empirical PSF targeted with a similar observing mode. 
Unfortunately, only a few bright stars exist in our targeted 
fields and we did not propose for separate PSF exposures. 
Comparison with a bright star in one field showed that 
the TinyTim model PSF reproduced the stellar PSF well 
(±5% in intensity; see Figure 2) out to a radius of about 
3". The final PSF-convolved template library consists of 
1000 model galaxies from which azimuthally-averaged pro¬ 
files have been extracted. Then a subset of 200 are chosen 
for each galaxy based on its measured ellipticity. 

We used a y^-minimization algorithm (AMOEBA; Press 
et al. 1992) to derive the sky value, which is an additive 
normalization, and a multiplicative normalization which 
gives the best fit to the data. Each image should be com¬ 
pletely sky-subtracted by the pedestal removal process (as 
the small fit sky values confirm); varying the sky value 
only lengthened the fitting process slightly while main¬ 
taining the greatest possible source-to-source consistency. 
We calculated the input errors for the procedure by 
measuring the output weight maps from DRIZZLE in the 
same elliptical isophotes as the galaxy. The best-fit model 
for each source was determined by comparing the resultant 
values for each of the 200 templates. The central «075 
diameter has been excluded from these fits to reduce the 
influence of any nuclear point sources and uncertainties 
associated with the PSF and centering. 


3.1.2. Two-Dimensional Fits 

Given a set of input parameters, the two-dimensional fit¬ 
ting algorithm generates model images “on-the-fly”. The 
freedom of the 2D approach allows us to fit for any combi¬ 
nation of galaxy centroid, position angle, ellipticity, effec¬ 
tive radius, surface-brightness, point-source contribution, 
and sky value. The method can be easily extended to any 
analytic profile. We primarily attempted to distinguish 
between de Vaucouleurs and exponential disk profiles, but 
also fit a Sersic law to determine whether deviations 
from the or disk models were prominent. We did not 
fit more complex profiles with additional free parameters 
because we felt it was not warranted by the signal-to-noise 
ratio and physical resolution of our data. 

The program extracts a thumbnail image centered on 
the galaxy from the larger, sky-subtracted mosaic. Ini¬ 
tially, this sub-image was 200 drizzled pixels (775) on a 
side, with the outer fit radius fixed to half this value (3775). 
Trials with an outer fit radius twice as large (775) did 
change the fits slightly for several galaxies. The average 
decrease in effective radius. Re, was « 10%. We chose one 
of these outer fit radii based on a simple by-eye compar¬ 
ison of the residuals for each source, rather than apply a 
strict test of xt values. 

The algorithm then refines the galaxy centroid to sub¬ 
pixel precision. An initial guess at the fit parameters is 
passed to a function that generates a model image which is 
then convolved with the TinyTim PSF (described above). 
The model is calculated on a sub-sampled pixel grid. Our 
sub-sampling by a factor of five allows the surface-brightness 
per pixel to be calculated to higher precision than at the 
original scale. This also allows centering to be accurate to 
less than a pixel. The model is then rebinned to the size 
of the input data thumbnail while conserving intensity. 

A xt is calculated for the model fit using the output 
weight map derived from the DRIZZLE task and a bi¬ 
nary mask generated to remove neighboring objects and 
to impose inner and outer fit radii. We masked neighbor¬ 
ing objects or high surface-brightness aligned light in the 
following galaxies: 3CR 65 (neighbor), 3CR 184 (neigh¬ 
bors), 3CR 210 (neighbors), 3CR 239 (small neighbors), 
3CR 266 (neighbors and aligned light), 3CR 267 (neigh¬ 
bors), 3CR 280 (neighbors and aligned light), and 3CR 368 
(cores of neighboring, largely subtracted stars and aligned 
light). The Levenberg-Marquardt minimization algorithm 
(Bevington & Robinson 1992) is used to generate the next 
guess for model parameters until the fractional difference 
between successive xt values is less than a small user- 
defined value. 

After testing several combinations of fit parameters (see 
§3.2) we also included a point source component in our 
fits and reduced the inner fit radius to zero. The fact that 
there is a degeneracy between the central cusp of a de Vau¬ 
couleurs profile and a nuclear point source, is discussed in 
more detail in §4.2. The simulations we performed con¬ 
vinced us of the importance of including this component. 
The galaxy centroid was allowed to move for all of the 2D 
fits. The point source centroid was always fixed to the 
galaxy centroid. In addition to analyzing our own data, 
we have used the two-dimensional method to fit 15 in¬ 
termediate redshift radio galaxies from the 3CR WFPC2 
snapshot survey (de Koff et al. 1996; Martel et al. 1999). 
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Fig. 3.— Percent deviations (output minus input) of derived parameters from input parameters as a function of input point source 
luminosity. The stars indicate fits which include a point source contribution in the models, the circles have only an -law component. 
These values were derived using the simulations fully described in § 3.2 


The results of the 2D fits are presented for the high- and 
intermediate-redshift RGs in Tables A and A and are dis¬ 
cussed in §4. 

3.2. Simulations and Error Analysis 

To better understand the systematic uncertainties in the 
hts, we undertook a thorough set of simulations. The full 
parameter space was explored by creating a large library of 
artificial galaxy images distinct from those used in the ID 
hts. Models were generated for a wide range of effective 
radius, nuclear point-source luminosity and signal-to-noise 
ratio, which is dehned per pixel at Re- We assumed a 
constant sky background, but the noise was modulated 
by the hnal weight-map produced by DRIZZLE (from a 
representative sub-image of a single target’s weight-map). 
These galaxies were then ht using the 2D method described 
above. Position angle and ellipticity were held constant at 
the same input value for all simulated sources. Table A 
shows the range of input values for the model galaxies. 

In Figure 3 we show the fractional deviation of the out¬ 
put ht values from the input values for radius, surface- 
brightness, and point-source strength versus point-source 
luminosity, which is shown as a percentage of the total 
galaxy luminosity. We can estimate the systematic errors 
on all the ht parameters by measuring the deviations of 
the output values compared to the input. Systematic er¬ 
rors on the primary ht parameters, effective radius, Re, 
and effective surface brightness, A, are less than 5% for 


all model galaxies with signal-to-noise ratio > 1 measured 
per pixel at Re- The median signal-to-noise ratio for our 
sources is w 3.5. Therefore, we expect the systematics on 
our derived parameters to be roughly in line with those 
simulated here. In Figure 3 stars represent hts which in¬ 
clude a point-source component, and circles represent hts 
which are purely de Vaucouleurs law models. The inclu¬ 
sion of a point-source in the models greatly reduces the 
systematic error on all three derived parameters: Re, le, 
and, of course, point-source luminosity. For point-source 
luminosities less than « 10% of the total galaxy light, the 
point-source hux itself is determined to w 5% of its in¬ 
put value. This is comparable to the intrinsic scatter for 
the galaxies with only a small point-source contribution. 
Therefore, the error in the nuclear luminosity determina¬ 
tion may introduce a spurious contribution to the total 
error budget when the nucleus is not very bright. 

Based on these simulations, we expect to underestimate 
the point-source contribution, and therefore underestimate 
the effective radius, because a more concentrated prohle 
accounts for the “missed” point-source light. The simu¬ 
lations can be used to quantify this effect as a function 
of point-source luminosity. For the derived point-source 
strengths of our targets, this is not a signihcant effect and 
no correction has been made. Therefore, we have chosen 
to perform hts to the data both with and without a point- 
source component, and determine which produces the best 
residuals, i.e., the lowest xt- 
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Effective Radius (pixels) 

Fig. 4.— contours representing one, two and three cr errors. The axes are linear intensity and radius, and the box is 15% of the best-fit 
value on a side. The best-fit value is marked with a cross. All of the fit galaxies share a similar x^ and error contour structure. The Icr error 
is 8% of the best-fit value in intensity and 5% of the best-fit value in radius. 


We measured the random errors using a technique. 
For each target galaxy we calculated x^ values for a fine 
linear grid of Re and le centered on the best fit values. All 
objects share essentially the same x^ structure. A repre¬ 
sentative contour plot is shown in Figure 4. The box is 
15% of the best fit values on each side, with radius the 
abscissa and linear intensity the ordinate. The contours 
represent the gaussian one, two, and three cr errors based 
on Ax^ values of 2.3, 6.17, and 11.8 respectively as cal¬ 
culated for the standard two degree-of-freedom table 
(Press et al. 1992). The median Icr error for all the sources 
is roughly 10% of the best-fit value in both radius and 
surface-brightness. 

The one- and two-dimensional fits agree reasonably well 
as illustrated in Figures 5 and 6. A linear regression anal¬ 
ysis calculates the scatter about the one-to-one relation to 
be 0"15 in radius and 0.3 mag in surface-brightness. Recall 
that these are completely independent methods of deter¬ 
mining the fit parameters. The agreement between the 
results of the two algorithms in addition to the results of 
the simulations gives us confidence that large, systematic 
errors in our derived galaxy parameters are unlikely. 



One—dimensional (orcsec) 

Fig. 5.— Effective radius derived via the 2D fitting method versus 
radius derived via the ID algorithm. The solid line is the one-to-one 
line. The two independent methods show good agreement with a 
scatter about the one-to-one line of 0(D5. 



Fig. 6 . — Effective surface brightness from the 2D fits versus the 
same from the ID fits. Again, the two independent fitting methods 
show good agreement. The scatter about the one-to-one line is 0.3 
magnitudes. 






Zirm, Dickinson & Dey 


11 


4. RESULTS AND DISCUSSION 

The NICMOS data presented here are unique in their 
combination of high spatial resolution and NIR wavelength 
coverage. Previous studies of z « 1 radio galaxies were 
limited either by contamination from the AGN at observed 
optical wavelengths, i.e., the UV alignment effect, or by 
ground-based seeing in the NIR. The exquisite resolution 
and NIR sensitivity provided by HST has allowed us to 
analyze these data in ways that were previously impossi¬ 
ble or unwarranted. In this section we discuss the results 
of our analysis, focusing first on the galaxy morphologies 
and then on their comparison to other galaxy samples. 
The key quantities used in the following discussion are 
presented in Tables A and A (and Table A for compari¬ 
son). Table A contains various magnitudes which are ei¬ 
ther model-independent or depend only weakly on the 2D 
fits, while Table A presents parameters derived directly 
from the fits which may depend on our choice of cosmo¬ 
logical model. 

4.1. Host Galaxy Morphology 

The NICMOS images (Figure 1) reveal round, symmet¬ 
ric host galaxies unseen in the UV WFPC2 data. Given 
the high angular-resolution of NICMOS, this unambigu¬ 
ously confirms the results of previous authors (e.g., Rigler 
et al. 1992; Dickinson et al. 1995; Best et al. 1998). The 
signal-to-noise ratio and spatial resolution of the data al¬ 
low a detailed investigation of the varied luminous compo¬ 
nents. In this section we present the analysis of the host 
galaxy morphologies. We studied the extracted profiles in 
detail to look for signs of both nuclear point sources and 
non-homology (i.e., the degree to which the profiles differ 
from one another). In combination with the comparably 
resolved UV data from WFPC2, the correspondence be¬ 
tween the AGN-related aligned light and the stellar host 
can also be determined. 

4.1.1. Stellar Component 

The one- and two-dimensional surface-brightness profile 
fits statistically prefer the models over the exponential 
disk profiles. To test this preference we have employed the 
statistical F-test to compare the resulting xl values. The 
disk model can be rejected in favor of the de Vaucouleurs 
profile at greater than 99.9% confidence for all objects. 
This suggests that the host galaxies are indeed ellipticals 
and provides circumstantial evidence that they are the pas¬ 
sively evolving galaxies responsible for the smooth K-z re¬ 
lation. The NIGMOS data presented here are the most di¬ 
rect detection of the hypothesized bulge component. The 
correlation between bulge and black-hole mass (Magorrian 
et al. 1998; Gebhardt et al. 2000; Ferrarese & Merritt 2000) 
also predicts that very luminous AGN should reside at the 
centers of large, bulge-dominated galaxies. We note, how¬ 
ever, that the de Vaucouleurs model is an oversimplifica¬ 
tion which ignores several details of the surface-brightness 
profiles. 

Figure 7 shows the smooth, simple surface-brightness 
profiles of the nine fitted galaxies. The solid line is the 
2D fit. Close inspection of the profiles reveals several 
features. Three galaxies, SCR 265, SCR 280 and SCR 65, 


show very steep central cusps, high central surface-brightness 
and the hint of the first diffraction ring of the NICMOS 
Camera 2 PSF (at r w 0"25). These clues suggest the 
presence of nuclear point sources or at least a very cuspy 
profile. This possibility we investigate more fully below 
and is included in the 2D fit for these S galaxies. SCR S68 
and SCR 210 show very flat outer profiles, seemingly more 
consistent with a power-law model rather than a de Vau¬ 
couleurs profile. Also, in both SCR 265 and SCR 280 the 
data exceed the profile at radii greater than « 077. 
This may be due to their extremely high central surface- 
brightness, because to compensate for the central surface- 
brightness, the law must be more compact and there¬ 
fore underpredicts the outer intensity. However, this ex¬ 
cess light at large radii may also be indicative of extended, 
cD-type halos (both 265 and 280 are believed to lie at the 
centers of moderately rich clusters), extended emission-line 
nebulae (for example, [Sill] A95SS in the case of 265) or 
the aligned light (particularly the “arc” to the west of 280, 
which was masked in the 2D fits but may still contribute 
measurable diffuse flux). 

Perhaps the more significant conclusion to be drawn 
from the differences between the galaxy profiles is that 
these sources may not form a homologous sample. Homol¬ 
ogy is relevant to the various scaling relations followed by 
early-type galaxies. For instance, the Fundamental Plane 
is tilted with respect to the predicted “virial plane” (e.g., 
Longo & Puddu 1997). One explanation for the tilt is non¬ 
homology among the galaxies. We have fit a Sersic profile 
to each source in order to test the assumption of homology 
for our sample. The Sersic profile (Sersic 1968) is defined: 

/(r) (1) 

where 6„ is calculated such that rg is the half-light radius. 
This more general profile has been shown to correlate with 
physical galaxy properties such as stellar velocity disper¬ 
sion (see Graham et al. 2001; Trujillo, Graham, & Caon 
2001). Sersic fits performed in the same manner as the 
and disk fits reveal a range of the index n within our 
sample: from 2.6 for 3CR 266 to 6.8 for 3CR 280 (see 
Table A). Due to the limited radial extent of the data, 
these results should be considered suggestive rather than 
conclusive, since the best-fit n value depends strongly on 
the behavior of the profile at large radii (particularly for 
n > 4). Deeper data will be able to distinguish between 
different n with greater confidence. Graham et al. (1996) 
have shown that brightest cluster galaxies often obey the 
flatter, n > 4 profile. It is interesting that several of our 
targets are also known to reside in moderately dense en¬ 
vironments, including 3CR 265 and 3CR 280, two of the 
galaxies exhibiting an n > 4 profile. Again, based on the 
statistical F-test, the Sersic fits do not significantly im¬ 
prove upon the r^!^ fits except in the case of 3CR 184 
where n = 2.6. The range of n within the sample may 
also partly explain the observed change in derived effec¬ 
tive radius when the outer fit radius is increased. There 
is a correlation between the derived Sersic index and the 
degree to which the change in outer fit radius affects the 
measured galaxy size. Table A (Col. 8) lists the Sersic 
indices for the galaxies. 







Radius (arcsec) Radius (arcsec) Radius (arcsec) 

Fig. 7. — Surface-brightness profiles with derived two-dimensional fits. The solid line is the best two-dimensional fit. The high central 
surface-brightness, diffraction rings and inclusion of a point-source in their best-fit models imply that SCR 65, 265 and 280 may harbor 
obscured nuclear sources. 


4.1.2. Host Galaxy Alignment 

Apart from the stellar host, the NICMOS data also con¬ 
tain other luminous components, some of which are di¬ 
rectly comparable to structures seen in the rest-frame UV 
images from WFPC2. The residuals shown in Figure 1 
are faint and often aligned with the radio axis, shown as 
lines in the middle panels of Figure 1. The direct cor¬ 
respondence between the NIR aligned residuals and the 
structures seen in the optical WFPC2 images suggests that 
much of the NIR alignment is simply the long-wavelength 
extension of the rest-frame UV alignment effect. The de¬ 
tailed photometric properties of the aligned light will be 
discussed in a future paper. 

Perhaps more interesting, however, is the observed align¬ 
ment between the stellar host and the radio axis. The 
luminosity-weighted position angle (PA) of the stellar host 
has been measured using the IRAF task ELLIPSE (see 
§3.1) with any bright aligned light excluded from the EL¬ 
LIPSE fitting (as described in § 3.1.1). The host galaxies 


are certainly less elongated than the UV continuum light, 
but surprisingly many remain within ~ 30 degrees of the 
radio axis (Figure 8). A KS test against a uniform dis¬ 
tribution shows that the APA distribution is inconsistent 
with a random distribution at the 99.9% confidence level. 
Of course, these measurements of PA depend on the ellip- 
ticity of the source. Very round galaxies, such as 3CR 65 
and 3CR 239, do not provide a meaningful PA. Best, Lon- 
gair, & Rottgering (1998) defined a metric, which they call 
“alignment strength”, Ug = e x (1 — This properly 

weights each measurement according to galaxy ellipticity 
(e) and APA. A positive value of as indicates alignment. 
The Og values for our sample are plotted as a histogram in 
Figure 9. The distribution has both positive mean equal 
to 0.07 and median equal to 0.06, indicating a tendency 
toward alignment. Two-thirds of the galaxies show posi¬ 
tive alignment strength, similar to the value of 70% found 
by Best et al. (1998). 

In order to determine whether or not the alignment is 
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APA (degrees) 

Fig. 8.— Difference between the position angles of the radio source and the host galaxy major axis, measured at optical rest-frame 
wavelengths. The top panel shows measurements from WFPC2 images of low-redshift radio galaxies, while the bottom shows NICMOS 
results for our high-redshift sample. The lower histogram is labeled with the target names. 


an intrinsic property of the stellar component or the result 
of contamination by the aligned component we performed 
simulations of randomly oriented galaxies with compara¬ 
ble aligned light. The aligned light component is treated 
using both a model and using the actual measured NIC¬ 
MOS residuals for 3CR 266. The results of these simu¬ 
lations suggest that both the ellipticity and PA measure¬ 
ments can be skewed by a modest amount of aligned light. 
Therefore, we believe that while these results are intrigu¬ 
ing, they could reflect some incomplete subtraction of the 
aligned light rather than a property of the stellar distri¬ 
bution. However, the same residual aligned light does not 
have much of an effect on the other results presented in 
this paper, contributing on average only 4% of the total 
light from the host galaxy. 

In addition, using our same software, we have measured 
ellipticities and derived sizes and luminosities for a sam¬ 
ple of intermediate redshift radio galaxies from the 3CR 
WFPC2 Snapshot Survey (de Koff et al. 1996; Martel et al. 
1999). The “snapshot” sources were chosen using sev¬ 
eral criteria: multiple exposures to allow cosmic-ray re¬ 
moval, extended galactic appearance (i.e., not quasars or 
“N” galaxies), no major dust features, and good signal-to- 
noise ratio (which eliminates many of the higher redshift 
sources). This selection resulted in 15 3CR galaxies be¬ 
ing chosen for fitting: 3CR 219, 234, 310, 314.1, 319, 332, 
346, 35, 381, 403, 436, 452, 63, 79, 93.1. The median 
redshift of the snapshot sample is 0.16 and the range is 
0.03 < z < 0.26. All the snapshot galaxies were observed 
with either the F675W or F702W. With a slightly negative 
mean Us, these sources do not show the same host galaxy 
alignment as the higher redshift 3CRs. We revisit these 
sources for more detailed comparison in Section 4.4. 


4.1.3. Dust 

The symmetric host galaxy is often centered on a lo¬ 
cal minimum of the UV light (e.g., 3CR 210, 3CR 368), 
suggesting that dust is playing an important role in the 
observed morphology of these sources (cf. Dickinson, Dey, 
& Spinrad 1995). In low-redshift radio galaxies, dust is 
quite common in the central region, often configured in ei¬ 
ther a disk or a dust lane (e.g., Martel et al. 1999; de Koff 
et al. 2000). In these local objects, the dust is known to be 
associated with the central engine. The physical resolution 
of our data is insufficient to make the same claim for our 
targets. However, the presence of molecular gas in some 
intermediate- to high-redshift radio galaxies (Papadopou- 
los et al. 2000) also implies that conditions suitable for 
dust exist. The existence of dust in these sources may also 
be inferred by the presence of scattered light, although 
electrons could instead be responsible. 

4.2. Nuclear Point Sources 

AGN unification theories propose a quasar nucleus is 
hidden within every radio galaxy. In the rest-frame ul¬ 
traviolet, the nucleus is heavily obscured along the line- 
of-sight by a dusty torus, but as one observes at longer 
(infrared) or much shorter (X-ray) wavelengths, the opti¬ 
cal depth toward the nuclear source decreases significantly. 
Therefore, our rest-frame optical NICMOS images may be 
able to detect nuclear point sources where shorter wave¬ 
length observations could not. To search for possible nu¬ 
clear point sources we added a delta-function at the center 
of the model galaxies which, upon convolution, becomes 
the PSF. This method is complicated by a degeneracy be¬ 
tween the central surface-brightness produced by the de 
Vaucouleurs model and the point-source. Also, the galaxy 
may not follow the r^!^ -law in the center, and instead 
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Fig. 9.— Alignment strength (as defined by Best et al. 1998) for the z ^ 1 and intermediate redshift radio host galaxies. Positive as 
indicates an alignment with the radio axis. The high-redshift distribution in the lower panel has both positive mean (vertical solid line) and 
median, indicating a tendency for the stellar host to align with the radio axis. The histogram is labeled with the target names for clarity. 
The lower redshift sample, in the upper panel, shows no such alignment, with a slightly negative mean (solid line) and median. 


have a more steeply rising, “cuspier” profile which would 
also mimic a central point-source. A possible resolution to 
this conflict is to deconvolve the data, an effort we have 
not undertaken at this time. The delta-function approach 
does allow us to separate the galaxies into those with and 
without point sources with fairly good accuracy as shown 
by the simulations (§3.2). 

As a corollary to this method, we can use results from lo¬ 
cal, well-resolved galaxies to potentially further constrain 
the nuclear flux contribution. Previous work on low-redshift 
AGN and “normal” ellipticals has broadly categorized surface- 
brightness prohles at small radii (Lauer et al. 1995). The 
profiles are well-described by a “broken” power-law. In 
general the slope becomes shallower within a “break” ra¬ 
dius of 100 - 300 pc and does not switch sign. If the 
same rule-of-thumb holds true for distant galaxies, then it 
is relatively easy to derive both an upper and lower limit 
to the contribution of a nuclear point source to the galaxy 
surface-brightness. At the redshifts of our sources, the 
PSF spans roughly 600 pc FWHM, so we are not able to 
detect the power-law break directly. However, a flux es¬ 
timate for a central point source component is greatly af¬ 
fected by the model behavior in the central portion of the 
PSF. If we Ht a single power-law and see a central excess in 
the residual image, we have a robust detection of a cuspy 
profile which implies the presence of a nuclear source. 
While this is a sound theoretical argument, in practice 
none of our sources exceed this extrapolated power-law. 
This does not rule out the presence of point sources, but 
suggests that their luminosities must be relatively low. Of 
the nine galaxies for which we Ht full 2D surface-brightness 
models, three of them show detections of nuclear point 
sources using the delta-function approach. None of the 
three point sources, in 3CR 65, 265, and 280, are brighter 


than w 5% of the total galaxy light. This is consistent 
with the null result of the power-law method. We have 
calculated robust la upper limits to the point-source con¬ 
tributions in the other six galaxies by gradually increasing 
the nuclear luminosity while allowing galaxy centroid. He, 
and le to vary, until increases to a certain value. These 
upper limits are included in column 6 of Table A. 

4.3. Galaxy Luminosities and Sizes 

Using the results of the 2D fits, the galaxy redshifts, and 
an assumed cosmology (§1), we derive absolute magnitudes 
in several bands and physical sizes from the model fits. We 
then compare these quantities to those measured in other 
galaxy samples. 

Table A presents the apparent magnitudes of the host 
galaxy calculated three ways: from the light within a 63.9 
kpc metric aperture for the data (Col. 3), from the light 
within a 63.9 kpc metric aperture for the model fit (Col. 4), 
and from the total light of the model fit (Col. 5). We also 
list the magnitudes of the nuclear point source component 
(Col. 6). These host magnitudes are largely free of contam¬ 
ination from AGN-related components. Indeed, we have fit 
and subtracted the nuclear point sources and the aligned 
light is faint. The degree to which these two components 
do affect the raw data (as measured by the aperture mag¬ 
nitudes) can be estimated in a mostly model-independent 
way. The point-source component is centered on the peak 
of the profile and the absolute upper limit to its contribu¬ 
tion is equal to twice the flux contained within a diameter 
equal to the PSF FWHM (because for a gaussian PSF half 
of the total flux is contained within the FWHM). Similarly, 
most of the aligned light can be measured in the residual 
image where both the host galaxy and point-source have 
been subtracted away. The magnitude of the residuals is 
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Fig. 10.— Size distributions for the samples described in the text; (a) Schade et al. (1999), (b) Jprgensen et al. (1995), (c) Graham et al. 
(1996), (d) Govoni et al. (2000) (FRIIs, solid histogram, FRIs, dotted histogram) and (e) and (f) the intermediate- and high-redshift 3CR 
galaxies from this work. The sizes of our sample galaxies seem to lie most closely to the low-redshift FRIIs in panel (d), however, the KS test 
statistically rejects this hypothesis. 


usually less than 4% of the total aperture galaxy magni¬ 
tude. The one outlier is SCR 368 which has « 14% of the 
total galactic flux in several bright aligned clumps. For all 
the galaxies the upper limit to the central flux is less than 
10% of the aperture flux. The total magnitudes quoted in 
column 5 of Table A are based on the model fits alone and 
are used throughout the following analysis. 

The NICMOS Alters, F165M and F160W, cover the rest- 
frame wavelengths near 6500 A in our sources. It is there¬ 
fore straightforward to correct our Flieo measurements to 
the rest-frame Gunn-r band using a model spectral energy 
distribution (SED) as described below. We have chosen 
rest-frame Gunn-r because of the wealth of comparison 
data in this band. The choice of SED makes little dif¬ 
ference for this small correction, but has a significant im¬ 
pact when we attempt to calculate magnitudes which are 
evolved to the present-day. Absolute un-evolved r mag¬ 
nitudes are presented in Tables A and A for the z w 1 
(NICMOS) and z « 0.16 (WFPC2 snapshot) SCR radio 
galaxies respectively. 

To calculate both rest-frame and passively evolved mag¬ 
nitudes we have adopted a model SED from the Bruzual- 
Charlot (BC) code (a 1999 version of Bruzual & Chariot 
1993; Bruzual 2000). The model SEDs in the BC sys¬ 
tem are f\ spectra, normalized to one solar luminosity. 
We have drawn our SEDs from the set of empirical single 
stellar population models with a Salpeter IMF and solar 
metallicity. To calculate the stellar age, we assume that 
our galaxies had their last major epoch of star-formation 
in an instantaneous burst at z/ = 2.5 and have evolved 
passively to their observed redshift. This choice of SED 
and stellar age implicitly determines galaxy colors, which 
we use to translate our NICMOS observations into other 
bandpasses, and mass-to-light ratios, which we use to cal¬ 


culate rough stellar masses for our sources. 

By integrating the model SED through the appro¬ 
priate NICMOS filter we normalized the spectrum to our 
observations. This normalized SED was then used to cal¬ 
culate the color correction to any other desired band. To 
correct to the Gunn-r band above we took the normalized 
SED, integrated it through the rest-frame r band filter 
and applied the AB to Gunn zeropoint offset. Then to 
calculate the rest-frame Gunn-r from the observed Ffieo 
we used: 

r = Hieo + 2.5 x logio((/i6o)/(A)) + 0.23 (2) 

where (/leo) and (fr) are the average fi, integrated over 
each filter and the 0.23 is the zeropoint offset going from 
AB to Gunn-r magnitudes (Fukugita, Shimasaku, & Ichikawa 
1995). Applying the same normalization to an evolved 
SED and calculating the corresponding color correction, 
we derive “evolved” magnitudes which we use to compare 
the z w 1 galaxies with local samples. 

We have also calculated rest-frame K magnitudes, in¬ 
cluded in column 10 of Table A. The more significant 
color correction from our observed wavelength to the rest- 
frame AT-band (also included in Table A) depends weakly 
on the assumed burst redshift via the amount of color evo¬ 
lution of the stellar population. If we change the redshift 
of the last major epoch of star-formation from 2.5 to 10, 
the mean absolute K magnitude brightens by 0.1 mag. 

Using the normalization and the gas mass fraction of the 
SED we also derive stellar masses. As shown in column 
11 of Table A these fall in the range of 10^^ Mq. These 
derived stellar masses are only weakly dependent on the 
assumed epoch of the last major episode of star-formation. 
We feel the choice of passive evolution is justified by both 
the K-z diagram and our later analysis of the size-surface- 
brightness relation (§4.4.3). 
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Fig. 11.— Galaxy luminosity distributions for the same samples as Figure 10. In panel (f) we have included both the observed luminosities 
of our sample galaxies (solid histogram) and those luminosities passively evolved to zero redshift (dotted histogram). The evolved magnitudes 
overlap the low and intermediate redshift radio galaxy distributions with relatively high statistical significance. 


The galaxy half-light radii are presented in Table A in 
both arcseconds (Col. 2) and kiloparsecs (Col. 3). For con¬ 
sistency with previously published studies, we have only 
used sizes derived using the -law profile rather than 
those from the Sersic fits. The galaxies are not perfectly 
round; the “radius” value we have quoted is the geometric 
mean of the semi-major and semi-minor axes. Radii for the 
2 w 1 3CR galaxies span the range from 3 to 17 kpc. We 
find the smallest galaxies also happen to be those with the 
largest point source contributions. We suspect that this is 
at least partly a coincidence, since our simulations show 
that the radius determination is accurate when fitting for 
the point source. 

4.4. Comparison With Other Samples 

To put our sample into the context of early-type galax¬ 
ies in general, we compare the z « 1 3CR galaxy sizes and 
luminosities to those drawn from previous studies. The ef¬ 
fective radii and host galaxy luminosities derived above are 
directly comparable to previous studies of r^!^ galaxies. 
For our comparison datasets, we have chosen the samples 
of Schade et al. (1999) (field early-type galaxies to z = 1), 
Jprgensen, Franx, & Kjaergaard (1995) (local cluster ellip¬ 
ticals), Graham et al. (1996) (local brightest cluster galax¬ 
ies), and Govoni et al. (2000) (local radio galaxies). We 
have also included the WFPC2 3GR snapshot survey data 
described above (§ 4.1.2) 

The Schade et al. (1999) sample was selected from the 
CFRS (Le Fevre et al. 1995) and LDSS (Colless et al. 1990) 
spectroscopic surveys based on two-dimensional surface- 
brightness fits to follow-up 77S'T/WFPC2 F814W imag¬ 
ing. Each galaxy included in their analysis was better fit 
by the law than either the exponential disk or a com¬ 
bination (disk-l-bulge) of the two models. The resulting 
catalog contains 48 elliptical galaxies and ranges in red¬ 


shift from 0.281 to 0.992. We cannot strictly compare the 
luminosities of this sample to our own because it contains 
galaxies over such a wide range of redshift. However, we 
feel that the more general comparison to such a sample has 
merit since the effective radii are still directly comparable 
to ours and the luminosities, while not strictly appropri¬ 
ate, should span the parameter space of all galaxies. The 
Jprgensen et al. (1995) study is a collection of data for 
211 early-type galaxies from 10 different nearby clusters. 
Photometry was performed on images obtained with the 
Danish 1.5m telescope at ESO. The targets were selected 
based on previous identification as cluster members. Mag¬ 
nitudes and sizes were derived using a one-dimensional 
fit to the photometric curve of growth. The Gra¬ 
ham et al. (1996) sample utilizes the photometry done in 
support of the Lauer & Postman (1994) study of peculiar 
velocities of Abell clusters. The fits to the BCGs were 
performed on isophotal profiles extracted from the ground- 
based images. The Govoni et al. (2000) sample is selected 
from two radio flux-limited samples (Wall & Peacock 1985; 
Ekers et al. 1989). Galaxy sizes and surface-brightnesses 
were derived via 2D image modeling of ground-based i?- 
band images of 79 low-redshift (z < 0.12) FRI and FRII 
radio galaxies. The authors include three components in 
their model fits, an r^!'^ profile, a nuclear point source 
(represented by the PSF) and an exponential disk. The 
three components were combined to obtain the optimal 
fit. For the size comparisons presented here we use the 
Govoni et al. (2000) empirically determined half-light ra¬ 
dius. For the luminosity comparison we have adopted their 
quoted total magnitudes, which are based on an extrapola¬ 
tion of the observed surface-brightness profiles rather than 
the model fits. However, for the purpose of constructing 
the Kormendy relation we have used solely the parameters 
from the de Vaucouleurs fit. It should be noted that the 
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Fig. 12.— Top panel: X-band luminosity function from Cowie et al. (1996) with arbitrary normalization. The histogram is the observed 
luminosity distribution of the z ^ 1 radio galaxies. The mean RG luminosity (solid line) falls at the very luminous end of the “normal” galaxy 
distribution (M* for the Cowie et al. (1996) and de Propris et al. (1999) samples are shown as the dotted and dashed lines respectively). 
Bottom panel: r'-band luminosity function from the Sloan Digital Sky Survey (Blanton et al. 2001) with arbitrary normalization. The 
histogram is now the passively evolved (to zero redshift) luminosity distribution of the z ^ 1 radio galaxies. Again, the evolved sources are 
brighter than M* (dotted line) suggesting that they are rare, massive sources. 


Govoni et al. RGs (particularly the FRIs) have total mag¬ 
nitudes (which include point-source and disk components) 
and sizes comparable to those of the BCGs. The galaxy 
luminosities and sizes of these samples have all been con¬ 
verted to our chosen cosmology. 

4.4.1. Sizes 

In Figure 10 we present the size distributions from our 
own sample along with the data from the literature de¬ 
scribed above. The NIGMOS galaxies have a wide size 
distribution spanning 3.7 to 17.1 kpc. To the eye, their 
sizes seem most similar to the FRII galaxies of the Govoni 
et al. (2000) sample (solid line in Fig. 10). We employ the 
KS test to statistically compare the z w 1 3GR galaxies to 
all the other samples. The local cluster, brightest cluster, 
and intermediate-redshift field galaxies are wholly incon¬ 
sistent with the NIGMOS 3GR size distribution at greater 
than 99.5% confidence. The most consistent score, at a 
statistically inconclusive 65%, is found for the local FRII 
galaxies as suggested by inspection. 

There is no exact counterpart population at low redshift. 
Lacking such a direct correspondence, we assert that it is 
unlikely that galaxy size decreases with time. Mergers will 
either leave the size unchanged, as in the case of accretion 
of a low-mass galaxy, or perhaps increase the size, in the 
case of equal mass mergers (e.g., Gapelato, de Carvalho, & 
Carlberg (1995)). If the merger funnels gas to the galaxy 
center, where it then forms stars, the merger may tem¬ 
porarily shrink the half-light radius. However, this effect 
will not dominate the observed light for long and is there¬ 
fore not likely to be a major contributor to observed galaxy 
properties. Under the scenario where their sizes either re¬ 
main the same or increase slightly, the z « 1 3CR galaxies 


match up with the local RGs or, if they grow, overlap the 
BCGs. 

The one outlier at small size (3CR 332 at « 0.3 kpc) in 
the WFPC2 snapshot sample is contaminated by a strong 
point source. This source has a doubly-peaked Ha line 
as shown by Halpern (1990), suggesting that we may be 
observing the nucleus/accretion-disk more directly than in 
most other radio galaxies. The derived nuclear luminos¬ 
ity in this case is much higher relative to the host galaxy 
than for the NIGMOS galaxies. If this galaxy is removed, 
the intermediate-redshift radio galaxy size distribution be¬ 
comes more consistent with the high-redshift one. 

4.4.2. Luminosities 

In Figure 11 we present the luminosity distributions for 
the NIGMOS sample along with the previously published 
samples described above. Like the size distributions, the 
z w 1 radio galaxy luminosities are roughly bracketed by 
the luminosities of the other samples. Again, we use the 
KS test to compare the distribution of our sample galaxies 
to those from the literature. In contrast to the sizes, how¬ 
ever, we perform this comparison for the observed values 
directly, but also for the observed values modified using 
an assumed evolutionary scenario. 

As determined by the KS test, the high-redshift radio 
galaxies are not drawn from the low-redshift cluster el¬ 
liptical or intermediate-redshift radio galaxy distributions 
at greater than 99.9% confidence. The exceptions are the 
BCG sample and the local FRI sample. The value for 
the BCG comparison is a statistically indeterminate 60%, 
meaning that the z « 1 3CR sample is inconsistent with 
being drawn from the BCG sample at this significance. 
The value for the low-redshift FRI sample is an also incon- 
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Fig. 13.— Kormendy relation in rest-frame Gunn-r. The dots are cluster ellipticals from Jprgensen, Franx, S>c Kjsergaard (1995). The 
pluses are the low-redshift galaxy sample of Govoni et al. (2000). The triangles are intermediate redshift (-^median = 0.16) 3GR radio galaxies 
from the WFPC2 snapshot survey (Martel et al. 1999), fit with our code. The crosses are the brightest cluster galaxies from Graham et al. 
(1996). The large squares are the NIGMOS 3CR data and are labeled by name. A representative average error-bar for the NIGMOS data is 
shown in the lower left panel. 


elusive, but slightly more consistent, 20%. It is interesting 
to note that while the sizes seem most consistent with the 
low-redshift FRIIs, the luminosities coincide with the lo¬ 
cal FRIs. In order to make a more direct comparison with 
local samples we have evolved the z « 1 SCR galaxies to 
z = 0 using the method outlined above (§ 4.3). 

When we do this, the luminosity distribution for the 
current sample shifts fainter by 1.5 - 2.0 magnitudes via 
passive evolution to zero redshift (dotted line in panel f 
of Figure 11). The amount of evolution was calculated 
using the normalized model SED described in § 4.3. We 
applied the same normalization to the model spectrum 
aged to zero redshift, and integrated it through the rest- 
frame Gunn-r filter. With this sizable correction the mag¬ 
nitudes are again bracketed by the low-redshift samples. 
Now the BCGs, cluster galaxies, and local RGs are re¬ 
jected at > 99% significance. While no sample is exactly 
comparable, the intermediate-z RGs are only rejected at 
the 34% significance level. If 3GR 332 is again thrown 
out on the basis of its discrepant size, the luminosity dis¬ 
tribution of the intermediate-redshift RGs is rejected at 


only the 25% level, making it the most consistent of the 
comparison samples. 

Another perspective can be taken by comparing the lu¬ 
minosities with much broader, more statistically robust 
samples. We can judge the relative rarity of the radio 
galaxies when compared to the more general luminosity 
functions (LFs) derived from large galaxy samples. At 
low redshift, the Sloan Digital Sky Survey provides a large 
number of galaxies for which the luminosity functions have 
been computed (Blanton et al. 2001). At higher redshift 
(0.6 < z < 1.0), we compared our 3GR galaxies to the 
field galaxy LF from Gowie et al. (1996) and the cluster 
galaxy LF from de Propris et al. (1999). The Gowie et al. 
(1996) spectroscopic targets were chosen from an infrared 
imaging survey covering 26.2 □' of the Hawaii Deep Field 
(areas SSA 13 and 22). We use their measurement of 
and a for the 0.6 < z < 1.0 bin as a comparison to our 
sample. The more recent study done by de Propris et al. 
(1999) investigates the AT-band LF in galaxy clusters to 
z « 1. In Figure 12 we show a histogram of our galax¬ 
ies along with the LFs from the Gowie et al. (1996), de 
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Propris et al. (1999), and (Blanton et al. 2001) samples. 
Using the assumed passive evolution model as above, we 
have calculated both the rest-frame iti-band (with no evo¬ 
lution) and r-band magnitudes (evolved to zero redshift) 
to compare with these two LFs. In both cases, the HzRGs 
are at the very luminous end of the LF, underscoring their 
relative rarity. 

4.4.3. Kormendy Relation 

Local spheroids (early-type galaxies and spiral bulges) 
define a continuous and low-scatter sequence in the size- 
surface-brightness plane. This two-dimensional projection 
of the three-dimensional Fundamental Plane (FP) provides 
a useful diagnostic of galaxy populations without velocity 
dispersion measurements. While the physical processes re¬ 
sponsible for the correlation are not entirely understood, a 
scaling relation of this sort is predicted qualitatively by the 
application of the virial theorem to these stellar systems 
(assuming a single mass-to-light ratio). Although this may 
only be illustrative, a galaxy which lies on the Kormendy 
relation is likely to be in dynamical equilibrium and may 
share a similar formation history with the other galaxies 
on the size-surface-brightness locus. 

Figure 13 shows our sample alongside the same low- 
redshift comparison galaxies, excluding the intermediate 
redshift field ellipticals from Schade et al. (1999). Of 
these samples, the local cluster, brightest cluster, and ra¬ 
dio galaxies are known to lie on the FP (Jprgensen, Franx, 
& Kjaergaard 1995; Oegerle & Hoessel 1991; Bettoni et al. 
2001). The 3CR snapshot galaxies form a tight sequence 
in the size-surface-brightness plane, but kinematic mea¬ 
surements are not yet available for their placement on the 
FP. A representative Icr error bar for the 2 ; « 1 3CR galax¬ 
ies, based on our error analysis in § 3.2, is shown in the 
lower left corner of the figure. 

The z « 1 3CR galaxies form a narrow sequence in the 
size-surface-brightness plot. This is not expected a priori 
for a set of sources selected via AGN rather than stel¬ 
lar properties (i.e., radio flux). To test the significance of 
this correlation defined by the high-redshift 3GR galax¬ 
ies, we have employed a simple rank correlation analysis. 
The small value of the two-sided significance indicates a 
strong correlation. The slope of the Kormendy relation is 
well-defined for low-redshift galaxies where evolutionary 
effects within the galaxy samples are small. For the NIG- 
MOS sources presented here this is not the case; several 
gigayears (= 3.3 for our cosmology) have passed between 
the observed epochs at z = 1.8 and z = 0.8. To allevi¬ 
ate this systematic bias we have evolved our sources to a 
common redshift, z = 0, following the procedure outlined 
above and assuming no evolution of galaxy sizes. 

Figure 14 again shows the same comparison galaxies ex¬ 
cept with the NIGMOS sources now passively evolved to 
zero redshift. The differential evolution between the sam¬ 
ples allows the z w 1 galaxies to fall roughly on the local 
relation. The slope of the HzRG relation is steeper (w 4.7) 
than that of the local galaxies (« 3.5). If we throw out 
the three point source galaxies (3GRs 65, 265, and 280) 
the slope is much shallower (w 3.6) and is more in line 
with the local value. Best et al. (1998) and McLure & 
Dunlop (2000) correspondingly derive slopes of 5 and 3.5 
for their similar sample of z « 0.8 RGs. It is interesting 


to note that our values range between these two slopes de¬ 
pendent upon the inclusion of the galaxies with a nuclear 
point-source. The slope is known to be very sensitive, 
particularly in these small samples, to observational selec¬ 
tion effects. To test more robustly whether the evolved 
sources do, in fact, align with the low-redshift relation, we 
have utilized the two-dimensional KS test (Peacock 1983; 
Fasano & Franceschini 1987). This test rejects the hy¬ 
pothesis that the z « 1 3CRs (evolved to zero redshift) 
are drawn from the same distribution as any of the local 
samples, at relatively high significance (81% for BCGs to 
99.9% for cluster members). This does not contradict the 
observed overlap between the high- and low-redshift sam¬ 
ples, but rather shows that they are not statistically identi¬ 
cal (e.g., as we have already shown, the size and luminosity 
distributions of the NIGMOS 3GR galaxies differ from that 
of any particular comparison sample). A more lenient test 
is to look at the surface-brightness distributions in the size 
range of the HzRGs. Using the one-dimensional KS test for 
these distributions, we find that the evolved 3GR galaxies 
are most similar to the low-z FRIs, followed by (in order of 
increasing rejection) the intermediate-redshift RGs (from 
the WFPC2 fitting), the low-z FRIIs, and the low-z cluster 
galaxies. In any case, the HzRG surface-brightnesses over¬ 
lap those for galaxies of similar size at low redshift when 
a simple model of passive luminosity evolution is applied. 

4.5. Radio Galaxy Evolution 

The NIGMOS data presented here are unique in their 
combination of high spatial resolution and near-infrared 
wavelength coverage. These capabilities have allowed us 
to study the properties of z « 1 radio host galaxies in some 
detail. The RG sizes, luminosities, and morphologies all 
have implications for the complete picture of radio, as well 
as perhaps “normal”, galaxy formation and evolution. 

A single, passively evolving population which formed at 
early times (z > 2.5) fits the observed K-z relation for 
HzRGs very well (e.g., Lilly & Longair 1984; Jarvis et al. 
2001). This result has defined the current radio host evo¬ 
lution paradigm: a massive, passively evolving stellar host 
dominates the rest-frame optical light and is responsible 
for the K-z relation, and a UV bright, but dynamically in¬ 
significant, component accounts for the alignment effect. 
Our data focus on the putative stellar component. Are 
the rest-frame optical morphologies, sizes, and luminosi¬ 
ties consistent with the passive evolution picture? 

The existence of a luminous, apparently stellar, con¬ 
tinuum component undetected in the rest-frame UV sup¬ 
ports the paradigm at its most basic level. Furthermore, 
the galaxy morphologies are well-described by the de Vau- 
couleurs profile, which at low redshift is associated with 
old, red elliptical galaxies. Ellipticals in clusters and the 
field seem to be passively evolving from early formation 
times (e.g., Stanford, Eisenhardt, & Dickinson 1998; Schade 
et al. 1999). In this context, the fact that the HzRGs are 
well-fit by the de Vaucouleurs profile suggests that they be¬ 
long to a similar, if not the same, class of passively evolving 
sources. 

The passive scenario also dictates that the radio hosts 
will undergo little or no size evolution. The HzRGs have 
similar effective radii to the low-redshift samples, again 
supporting the paradigm: the sizes of the z « 1 RGs are 
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Fig. 14.— Same as Figure 13 except the NICMOS galaxies have been passively evolved to zero redshift using the stellar population model 
(^Zf = 2.5) described in the text. The z ^ 1 3CR galaxies evolve roughly onto the local relation, but are not statistically consistent with the 
distributions of the comparison samples. 


most consistent with the z « 0 FRII radio hosts. We may, related light boosting the derived values. Based on the fits 

therefore, identify the HzRGs as the progenitors of local and residual images, we believe that the AGN contribution 

radio galaxies. Radio galaxies would then be drawn from cannot account for a majority of the luminosity difference, 

a single, passively evolving population of sources formed We therefore conclude that the z w 1 RGs are inherently 

at a common high redshift. The relative rarity of RGs at more luminous. If we assume they dim passively, then 

any given epoch is explained by noting that the radio emis- they fall in the range of local galaxy luminosities. Impos- 

sion in FRII sources is likely to be a transient phenomenon ing the galaxy radius restriction as above, the high-z 3GR 

occurring on timescales of roughly 10^ years (e.g., Alexan- galaxies are fainter than the BGGs and low-redshift RGs, 

der & Leahy 1987; Mack et al. 1998; Kaiser & Alexander and are brighter than the local cluster galaxies. In partial 

1999). To follow this argument further, we now discuss support of the paradigm, they line up most closely with 

the luminosities of the HzRGs. the intermediate-redshift radio hosts. 

With no evolution, the z w 1 radio galaxies are much In the size-surface-brightness plane the HzRGs evolve 

more luminous than all but the BGGs and BGG-like low- onto the local relation, although distinctly distributed from 

redshift FRI RGs. This is not a completely fair comparison any of the comparison samples. Again the HzRG distri- 

since BGGs are perhaps the product of a “special” forma- bution is most consistent with the local RGs. This cor- 

tion history (e.g., West 1994; Aragon-Salamanca, Baugh, respondence, like the evolved luminosities above, depends 

& Kauffmann 1998; Burke, Gollins, & Mann 2000) and on an assumed star-formation history. In order to fade 

are known to have distinct size and luminosity distribu- these sources onto the local Kormendy relation and into 

tions. If instead we restrict the size range of the compari- the luminosity range of local galaxies using the single-burst 

son samples to the same set of values spanned by the high- model, the most recent star-formation cannot occur too 

redshift 3CRs, then the z « 1 sample outshines all other early, roughly z/ w 2 — 2.5. However, most of our galaxies 

galaxies including the BGGs. This discrepancy could be would evolve sufficiently with z/ < 10. It is the highest 

due to inherently more luminous galaxies or to the AGN- redshift source in our fit sample, 3CR 239 at z = 1.78, 
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that drives to lower redshift the choice for the last ma¬ 
jor star-formation (which is not necessarily the same as 
the last major merging epoch). For a galaxy at z = 1.8, 
changing the “formation” redshift from 2.5 to 5 decreases 
the amount of fading to z = 0 by about 0.8 magnitudes. 
For the other, lower redshift, galaxies in our sample the 
difference is smaller, about 0.3 magnitudes. Such a low 
“formation” redshift to accommodate SCR 239 is prob¬ 
lematic because the RG K-z relation is known to extend 
to z > 5 with low scatter. In addition, SCR 239 exhibits 
several discrete clumps of emission surrounding the central 
concentration, suggesting possible residue of formation or 
in-falling satellites. It is also interesting to note that the 
highest redshift object in the entire sample, SCR 256 at 
z = 1.82, shows no centrally peaked emission and is un- 
derluminous at K compared to the K-z relation at that 
redshift (Chokshi & Eisenhardt 1991; Simpson et al. 1999). 
While two data points do not constitute a robust sample, 
these galaxies suggest that the epoch for radio galaxy for¬ 
mation may extend in some cases to z < 2. At z > 2 it 
has been shown that a majority of RGs, while falling on 
the AT-band Hubble diagram, show few signs of smooth, 
round host galaxies (van Breugel et al. 1998; Pentericci 
et al. 2001). So, either we maintain that passive evolution 
is consistent with our data and 3CR 239 is an outlier (al¬ 
though, unlike 3CR 256, it does not deviate significantly 
from the mean K-z relation), or we must imagine an al¬ 
ternate explanation for the observed Hubble diagram. 

If the distant radio hosts do not evolve into the local 
radio-loud galaxies, how is the K-z relation explained? An 
alternative posited by Best et al. (1998) is that low- and 
high-redshift RGs are not linked by evolution, but rather 
that RGs are only created in host galaxies above a given 
large galaxy mass. Their main supporting evidence for this 
view was that the cluster environments of RGs seemed to 
become less rich at lower redshift, something that cannot 
happen. The inclusion of a galaxy mass cutoff makes sense 
in terms of the relationship which has been discovered at 
low-redshift between black-hole mass and bulge velocity 
dispersion or luminosity (Magorrian et al. 1998; Gebhardt 
et al. 2000; Ferrarese & Merritt 2000). In this scenario, 
the natural relation between radio luminosity and redshift 
in a flux-limited sample populates the Hubble diagram 
in a smooth, continuous way via the correlation of radio 
power with black hole and bulge mass. While our data 
cannot address the environment issue, we do not see any 
comparable “impossible” transformation (e.g., low-z RGs 
being smaller than HzRGs or fainter than passive evolu¬ 
tion would predict) that precludes radio galaxies at high-z 
from being drawn from the same population as low-z RGs. 
Given our data, the simplest explanation of the K-z dia¬ 
gram, passive evolution, remains tenable. 

The fact that RGs form a highly-correlated group and 
follow largely the same analytical profiles (with possibly 
interesting differences) suggests that they may form a ho¬ 
mogeneous sample. Similar behavior among local galaxy 
samples has been used to argue that the constituent galax¬ 
ies have followed similar evolutionary paths. While our 
sample is too small to make such grandiose statements, the 
size-surface-brightness data certainly do not preclude the 
possibility that this radio-selected sample constitutes a ho¬ 
mogeneous subset of the entire galaxy population. This is 


seemingly not true for the higher redshift objects, 3GR 239 
and 256 in our sample and the z > 2 galaxies of van 
Breugel et al. (1998) and Pentericci et al. (2001), which 
begin to show interesting diversity in their properties. 

5. CONCLUSIONS AND FUTURE 

Historically, there have been many attempts to unify 
the disparate qualities of high-redshift radio galaxies into 
a full picture of their formation and evolution. KST and 
NIGMOS now play a crucial role in this continuing story. 
The present observations have targeted 11 z w 1 3GR radio 
sources for detailed study. Round, smooth galaxies appear 
in the NIR and outshine the complex morphologies seen 
in the rest-frame UV. The NIGMOS data have allowed 
us to analyze these sources in ways that were previously 
impossible or at least ill-advised. 

The salient results of our investigation are briefly sum¬ 
marized: The majority of our sample are well-fit by the 
de Vaucouleurs profile, with low-level residuals that qual¬ 
itatively match the rest-frame UV morphologies observed 
with WFPG2. The derived surface-brightnesses and sizes 
match well with bright local ellipticals, but are both fainter 
and smaller than the brightest galaxies found in rich clus¬ 
ters at low-redshift. Apart from one fit source (3GR 239), 
passive evolution from an early formation epoch (z > 5) 
is sufficient to fade the sources to match the local sam¬ 
ples. The slope of the Kormendy relation for our sample 
is steeper than that found for local galaxies, 4.7 compared 
to 3.5, and may be consistent with a constant luminosity 
line. If a lower “formation” redshift is chosen, z = 2.5, 
3GR 239 can simultaneously be made to fade into the 
range of low-redshift galaxies. However, comparison to 
luminosity functions at both low and high redshift illus¬ 
trates the relative rarity of galaxies at these luminosities, 
and suggests that HzRGs are quite massive. The z w 1 
3GR host galaxies may show a tendency to align their ma¬ 
jor axes with the axis defined by the double radio lobes, 
but confirmation of this result awaits detailed multi-band 
analysis. 

In short, the z « 1 radio galaxies seem to be rare, mas¬ 
sive objects which are dynamically relaxed. As such they 
are likely the progenitors of local giant elliptical galax¬ 
ies. These sources are ideal for follow-up NIR spectroscopy 
on large telescopes. They are relatively bright at H and 
should have measurable velocity dispersions at relatively 
low spectral resolution. In addition, similar studies of 
higher redshift RGs will test whether the evolutionary 
trends seen here between z = 0 and z w I continue to 
even earlier epochs. 

APPENDIX A 

NOTES ON INDIVIDUAL OBJEGTS 

SCR 265 (z = 0.811): This galaxy lies in a moderately 
rich cluster of galaxies (Deltorn 1998) and is surrounded by 
spectacularly extended UV emission seen in the WFPG2 
images. The extended UV light is strongly polarized (Jan- 
nuzi & Elston 1991) and shows extended, broad Mgll emis¬ 
sion (Dey & Spinrad 1996), implying that it is dominated 
by scattered AGN light. It is, however, significantly mis¬ 
aligned from the radio source axis. The NIGMOS image is 
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wholly dominated by a round, symmetric giant elliptical 
host galaxy. The central surface-brightness is remarkably 
high, and the model subtraction shows extended, irregular 
residuals, harboring an obscured point source. The point 
source fit hnds a central nucleus of rript^ieo = 22.84 in the 
F160W hlter. 

SCR 184 (z = 0.994): Another cluster galaxy (Deltorn 
et al. 1997), SCR 184 is elongated, bimodal, and aligned 
in the rest-frame UV, but has a dominant gE host galaxy 
in the NICMOS image. The isophotal elongation of the 
NICMOS image along the radio source axis is most likely 
due to the low-level, aligned light seen in the residuals. 
The WFPC2 image here is from Deltorn et al. (1997). 

SCR 280 (z = 0.996): The complex aligned structure 
visible in the WFPC2 images of Best et al. (1997) largely 
disappears in the NICMOS data, although the highest sur¬ 
face brightness components remain barely visible. The ra¬ 
dio galaxy host appears to be a smooth, symmetric giant 
elliptical, with a central point source (rript^ieo = 23.03). 

SCR S68 (z = 1.131): The rest-frame UV light from 
this object is extended over 7" in length along the ra¬ 
dio source axis and, surprisingly, is unpolarized. There is 
a galactic M-star superimposed on the galaxy (Hammer, 
Proust, & Fevre 1991), with another, fainter star at the 
north end visible in the NICMOS image. Both stars were 
subtracted before performing the surface-brightness fits. 
The host galaxy nucleus seen by NICMOS lies at a local 
minimum in the WFPC2 images; here and in other radio 
galaxies it is likely that a combination of dust-obscuration 
and fc-correction are responsible for rendering the stellar 
nucleus completely invisible at optical wavelengths. The 
near-IR central surface brightness of 3CR 368 is low and its 
surface brightness proHle is unusually extended and “flat”. 
It may be that dust obscuration remains significant in this 
object even at 0.8/xm in the rest frame. It is also possi¬ 
ble, based on the Sersic law fits, that the host galaxy of 
3CR 368 has a more BCG type profile (n > 4). 

SCR 267 (z = 1.140): The WFPC2 images of 3CR 267 
show diffuse emission, several bright knots, and a broad 
“arc” all distributed generally along the radio source axis. 
In the NICMOS image, the three brightest knots appear 
to be three distinct galaxies: the gE radio source host 
plus two fainter companions, one evidently a disk galaxy. 
3CR 267 is also known to inhabit a rich cluster. The fit 
residuals are very small in this case, without including any 
point source component. 

SCR 210 (z = 1.169): The complex, multicomponent 
structure seen in the WFPC2 data simplifies to a domi¬ 
nant, if not completely smooth, central host galaxy in the 
NICMOS image, together with another red companion ob¬ 
ject. Careful comparison of the WFPC2 and NICMOS 
data strongly suggests the presence of a thick dust lane 
obscuring the central regions in the rest-frame UV. Ex¬ 
tinction may be responsible for the central distortions still 
seen at the 0.8/xm rest frame in the NICMOS data. The 
3CR 210 host galaxy seems to have two distinct peaks, 
which are less distinct in the NICMOS data than in the 
rest-frame UV. This source also lies near the center of a 
moderately rich cluster (J-M. Deltorn, private communi¬ 
cation) . 

SCR 65 (z = 1.176): This galaxy was long believed 
to be a giant elliptical on the basis of ground-based near- 


IR images (cf. Rigler & Lilly 1994) and spectroscopy of 
the 4000A region (Stockton, Kellogg, & Ridgway 1995), 
although Lacy et al. (1995) have argued that it harbors an 
optically-obscured AGN point source which is detectable 
in the near-IR. Both points of view may be true: viewed 
with NIGMOS the galaxy appears to be a round gE, but 
with extremely high central surface brightness. The fits 
find a central point source with rriptpes = 22.75. 

SCR 266 (z = 1.275): One of the most elongated of 
all 3CR galaxies in the rest-frame UV, the alignment of 
SCR 266 persists into the near-IR, although the isophotes 
are generally rounder. The near-IR morphology is striking; 
cleanly separating into a relaxed, elliptical host and a thin, 
aligned component nearly identical to the WFPC2 image. 

SCR 470 (z = 1.653): This galaxy is extremely faint 
at optical wavelengths and is barely visible in the WFPC2 
images. Indeed it is not entirely clear which of the faint 
blobs in the WFPC2 data is the actual host: Best et al. 
(1997) place the radio nucleus between the two. Eisen- 
hardt et al. (1989) noted that the near-IR structure of 
SCR 470 is apparently oriented perpendicular to the radio 
source axis. In the NICMOS images it appears that this 
is due to the (chance?) alignment of three independent 
objects distributed roughly in a line. 

SCR 2S9 (z = 1.781): This galaxy is compact and 
bright in both the optical WFPC2 and NIR NICMOS 
images. The NICMOS image is dominated by a single, 
round, centrally concentrated galaxy, but its outer enve¬ 
lope breaks up into several “satellite” lumps. The image is 
suggestive of a giant but not fully mature elliptical galaxy 
at z w 1.78, perhaps seen in the process of assembling 
through accretion. While the proHle is well fit by an 
-law, the residuals again show that extinction may play a 
large role at rest wavelengths of 0.6/im. 

SCR 256 (z = 1.819): The most distant galaxy in our 
sample, 3CR 256 is extremely bright in the rest-frame UV 
and consists of two dominant components, strongly polar¬ 
ized, and closely aligned with the radio source (Jannuzi 
et al. 1995; Dey et al. 1996). It is elongated/aligned in 
the near-infrared, and is very subluminous, lying well be¬ 
low the mean K-z relation at that redshift. On this basis, 
Eisenhardt & Dickinson (1992) suggested that 3CR 256 
may be one of the few true examples of a protogalaxy 
among the 3CR radio sources (see (Simpson et al. 1999) 
for a detailed analysis of the spectral energy distribution). 
In our NICMOS F165M images the galaxy is diffuse, elon¬ 
gated and very faint: it is the one complete exception to 
the trend toward compactness and symmetry among the 
NICMOS 3CR observations presented here. One peak of 
the infrared light lies between the bright UV lobes; there is 
also a faint, compact red object barely visible in the optical 
data. An F160W NICMOS image (not shown here) which 
includes [OIII] line emission has a completely different ap¬ 
pearance, dominated by two bright lobes corresponding to 
those seen in the optical continuum and Lyman a. 
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Massive Ellipticals at z w 1 


Table 1 

Observation Log 


SCR 

RA®- 

Dec^ 

RA^ 

Dec^ 

Redshift 

Filter 

Rest-frame 

Exposure 

Dithers 

Estimated Line 


(J2000) 

(J2000) 

(J2000) 

(J2000) 



A Range (A) 

(s) 


Contamination*" 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 


265 

Ilh45m31.87s 

-|-31d33m43.4s 

Ilh45m29.0s 

-|-31d33m47s 

0.811 

F160W 

7805 

- 9997 

2559 

4 

< 0.06 

mag ([SIII]A9533) 

184 

07h39m24.47s 

-|-70d23ml0.9s 

07h39m24.2s 

-|-70d23m27s 

0.994 

F160W 

7089 

- 9080 

5134 

6 

< 0.01 

mag 

280 

12h56m57.49s 

-|-47d20m20.2s 

12h56m57.7s 

-|-47d20m21s 

0.996 

F160W 

7082 

- 9071 

5134 

6 

< 0.01 

mag 

368 

18h05m06.36s 

-|-Ild01m32.6s 

18h05m06.4s 

-t-lld01m33s 

1.131 

F165M 

7293 

- 8218 

10268 

12 

< 0.01 

mag 

267 

Ilh49m56.53s 

-|-12d47ml9.0s 

Ilh49m56.6s 

-f-12d47ml9s 

1.140 

F165M 

7263 

- 8183 

10268 

12 

< 0.01 

mag 

210 

08h58m09.96s 

-|-27d50m51.6s 

08h58ml0.0s 

-|-27d50m54s 

1.169 

F165M 

7166 

- 8074 

10268 

12 

< 0.01 

mag 

65 

02h23m43.19s 

-|-40d00m52.5s 

02h23m43.4s 

-H40d00m53s 

1.176 

F165M 

7142 

- 8048 

7701 

9 

< 0.01 

mag 

266 

Ilh45m43.37s 

-|-49d46m08.2s 

Ilh45m43.4s 

-l-49d46m08s 

1.275 

F165M 

6832 

- 7698 

8727 

13 

< 0.01 

mag 

470 

23h58m35.34s 

+44d04m38.9s 

23h58m36.0s*^ 

+ 44d04m45s*^ 

1.653 

F160W 

5328 

- 6824 

7701 

9 

w 0.20 

mag (Ha) 

239 

10hllm45.42s 

-|-46d28ml9.7s 

10hllm45.3s 

-l-46d28m20s 

1.781 

F160W 

5083 

- 6510 

7701 

9 

w 0.07 

mag ([OIII]AX4959, 5007) 

256 

Ilh20m43.05s 

-|-23d27m55.0s 

Ilh20m43.0s 

-l-23d27m54s 

1.819 

F165M 

5513 

- 6212 

7701 

9 

< 0.01 

mag 







F160W 

5014 

- 6422 

5134 

6 




®'Radio source coordinates (J2000) from Douglas et al. (1996) 

^NICMOS galaxy coordinates (J2000) as measured from image WCS, subject to absolute errors of ss 
*^Calculated assuming the composite spectrum of McCarthy (1993) 

‘^Assuming the middle source in the ‘trio’ is the radio galaxy. See the cutout in Figure 1 
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Table 2 

Galaxy Photometry'^ 


Object 

Filter 

«160'" ^ 
(Aperture Data) 

Hi60° 

(Aperture Model) 

(Total Galaxy) 

(Pt. Src.) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


SCR 265 

F160W 

18.67 

18.70 

18.74 

22.84 

SCR 184 

F160W 

19.27 

19.28 

19.04 

>30.03 

SCR 280 

F160W 

19.53 

19.51 

19.49 

23.03 

SCR 368 

F165M 

20.10 

20.18 

19.75 

>28.31 

SCR 267 

F165M 

20.19 

20.22 

20.07 

>25.12 

SCR 210 

F165M 

19.71 

19.68 

19.12 

>28.09 

SCR 65 

F165M 

19.26 

19.43 

19.35 

22.75 

SCR 266 

F165M 

20.09 

20.16 

19.92 

>29.07 

SCR 470 

F160W 

21.15^ 




SCR 239 

F160W 

19.88 

19.87 

19.68 

>27.19 

SCR 256 

F165M 

20.76 





®'A11 magnitudes are presented in the AB system 

^Measured in a 63.9 kpc diameter aperture, calculated in our adopted cosmology 

*^Total aperture light of the best fit de Vaucouleurs model and point source (where applicable) 

*^Total light in the best fit de Vaucouleurs model (no point source component) 

®Point source magnitudes. Upper limits are derived from the maximal allowed (within 1 <t errors) 
point source component, and are shown for galaxies which were better fit by the pure galaxy model 

^Assuming the middle source in the ’trio’ is the radio galaxy. See the cutout in Figure 1 
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Table 3 

Derived NICMOS Galaxy Parameters'^ 


Object 

(1) 

Outer Fit 
Radius {") 

(2) 

fie” 

(") 

(3) 

fie° 

(kpc) 

(4) 

(De” 

(mag/D") 

(5) 

Ellip.® 

(6) 

PA^ 

(°E of N) 

(7) 

Sersic Index 

(") 

(8) 

L. g 

ncovr 

(mag) 

(9) 

Me” 

(mag) 

(10) 

Mk‘ 

(mag) 

(11) 

Stellar Mass-^ 
(10” Mo) 
(12) 

3CR 265 

7.5 

0.46 

3.77 

17.87 

0.157 

171 


5.5 

0.52 

-23.58 

-24.57 

5.83 

3CR 184 

3.8 

1.10 

9.46 

19.59 

0.227 

82 


2.6 

0.37 

-23.86 

-24.85 

6.33 

3CR 280 

7.5 

0.55 

4.75 

18.53 

0.173 

91 


6.8 

0.37 

-23.42 

-24.41 

4.23 

3CR 368 

7.5 

1.54 

13.62 

20.75 

0.157 

171 


5.7 

0.32 

-23.49 

-24.45 

4.01 

3CR 267 

3.8 

0.77 

6.83 

19.53 

0.155 

49 


4.2 

0.30 

-23.21 

-24.11 

2.92 

3CR 210 

7.5 

1.93 

17.14 

20.50 

0.367 

10 


2.8 

0.27 

-24.24 

-25.14 

7.54 

3CR 65 

3.8 

0.68 

6.04 

18.45 

0.017 

8 


4.4 

0.26 

-24.02 

-24.92 

6.18 

3CR 266 

3.8 

1.09 

9.82 

19.83 

0.324 

5 


2.6 

0.19 

-23.70 

-24.57 

4.20 

3CR 239 

3.8 

0.90 

8.19 

18.16 

0.182 

180 


3.9 

-0.18 

-24.98 

-25.56 

5.90 


®'A11 magnitudes are presented in the AB system 
^Effective radius in arcseconds 
*^Effective radius in kiloparsecs 

Average surface-brightness within effective radius in rest-frame Gunn-r magnitudes per square arcsecond 
®Ellipticity 
^Position angle 

^fc-correction to rest-frame Gunn-r magnitude using our model SED described in § 4.3, 2.51og2^o((/i6o)/{/r)) 
Absolute (total) Gunn-r magnitudes 
^Absolute (total) rest-frame K magnitudes 

•^Stellar Mass as derived from the SED normalization to the NICMOS data 
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Table 4 

Derived WFPC2 Galaxy Parameters'^ 


Object 

(1) 

Redshift 

(2) 

Filter 

(3) 

Re*’ 

(") 

(4) 

Re° 

(kpc) 

(5) 

(De-* 

(mag/n") 

(6) 

Ellip.® 

(7) 

PA*’ 

(°E of N) 
(8) 

^corr 

(mag) 

(9) 

Mr** 

(mag) 

(10) 

3CR 35 

0.03 

F702W 

5.87 

3.85 

21.07 

0.281 

112 

0.08 

-20.42 

3CR 310 

0.05 

F702W 

10.57 

11.86 

21.58 

0.284 

81 

0.04 

-22.36 

3CR 403 

0.06 

F702W 

8.35 

10.27 

20.86 

0.320 

34 

0.02 

-22.77 

3CR 452 

0.08 

F702W 

3.54 

5.83 

20.23 

0.263 

103 

-0.02 

-22.17 

3CR 314 

0.12 

F702W 

2.70 

6.27 

20.68 

0.154 

76 

-0.10 

-21.87 

3CR 332 

0.15 

F702W 

0.11 

0.32 

14.85 

0.053 

54 

-0.16 

-21.24 

3CR 381 

0.16 

F702W 

1.50 

4.48 

19.54 

0.098 

53 

-0.18 

-22.28 

3CR 346 

0.16 

F702W 

0.44 

1.32 

17.99 

0.067 

113 

-0.19 

-21.18 

3CR 219 

0.17 

F702W 

0.79 

2.51 

18.41 

0.113 

110 

-0.21 

-22.16 

3CR 63 

0.17 

F702W 

0.71 

2.28 

18.43 

0.215 

84 

-0.21 

-21.92 

3CR 234 

0.18 

F675W 

0.74 

2.45 

18.32 

0.096 

109 

-0.30 

-22.20 

3CR 319 

0.19 

F702W 

0.58 

2.01 

18.70 

0.148 

152 

-0.25 

-21.38 

3CR 436 

0.21 

F675W 

4.39 

16.49 

21.41 

0.218 

11 

-0.36 

-23.25 

3CR 93 

0.24 

F675W 

1.37 

5.67 

20.28 

0.259 

130 

-0.43 

-22.06 

3CR 79 

0.26 

F675W 

1.59 

6.81 

19.40 

0.074 

20 

-0.45 

-23.34 


®'A11 magnitudes are presented in the AB system 
^Effective radius in arcseconds 
*^Effective radius in kiloparsecs 

^Average surface-brightness within effective radius in rest-frame Gunn-r magnitudes per square arcsecond 

®Ellipticity 

^Position angle 

^fc-correction to rest-frame Gunn-r magnitude using our model SED described in § 4.3, 2.51og2^Q({/i6o)/(/r)) 
^Absolute (total) Gunn-r magnitudes 
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Table 5 

Simulated Galaxy Parameters 


Parameter 

Values 

Re (pixels)®" 

3, 5, 10, 15, 30, 60 

Iq (cts/s/pixel) 

100 

Point Source Strength 


(% ot total galaxy light) 

0, 1, 5, 10, 20 

Signal-to-Noise 


(per pixel at Re) 

0.2, 1, 5, 15 


^These sizes correspond to effective radii of 
roughly 1-20 kpc at 2 = 1 in our assumed cos¬ 
mology, — (0.3, 0.7) and Ho — 65 



